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Introduction 

There were 5 sections of contributed papers on 
muons and neutrinos with a total of 107 papers 
distributed as shown in Table 1. The most active 
category is HE 2.3, astrophysical neutrinos. Par- 
ticularly in this area, there were also many papers 
presented in OG 2.5 sessions, on high-energy as- 
trophysical neutrinos. I include discussion of these 
topics to the extent necessary to present a coherent 
overview of the field as of mid-year 2007. 
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Table 1 : Papers on muons and neutrinos 







Session 


Topic 


# 


HE 2.1 


Muon experiments 


17 


HE 2.2 


Observations of solar 






& atmospheric v 


16 


HE 2.3 


Observations of 






astrophysical v 


37 


HE 2.4 


Theory and simulations 


19 


HE 2.5 


New experiments 






& instrumentation 


18 



Atmospheric muons 

Muons are the gold standard of cosmic -ray physics 
because they are well-measured and their physi- 
cal origin in the atmospheric cosmic -ray cascade 
is well-understood. A summary of measurements 



Figure 1: Muon charge-ratio from Ref. [4]. The 
line "vri^" model corresponds to a fit to the ratio of 
/Lt+//Lt^ fromEq. 1. 



is included in this conference in Ref. [1]. Muons 
are penetrating and relatively abundant in all ter- 
restrial particle detectors. They are therefore a po- 
tential source of background, and at the same time 
they are useful for detector calibration. One use 
of cosmic -ray muons is as a survey tool, some- 
times called muon tomography. A classic example 
is the survey of the Second Pyramid of Giza and 
the search for hidden chambers [2]. The status of 
a similar investigation of the Pyramid of the Sun at 
Teotihuacan was presented at this conference [3]. 
The detector is integrated and ready for installa- 
tion during 2008 in the tunnel that goes under the 
pyramid. 
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An important new result presented at this confer- 
ence is the measurement of the muon charge ra- 
tio in the far detector of MINOS [4, 6, 5]. Muons 
that reach the detector at its depth of 2070 meters- 
water-equivalent (m.w.e.) have energies at the sur- 
face in the range of 1 -7 Te V, depending on zenith 
angle. Fig. 1 from Ref. [4] shows muon charge 
ratio increasing in the energy range from i?^ < 
100 GeV to E/j, > 1 TeV. The potential signif- 
icance of this result can be understood from the 
energy-dependence of Eq. 1. 

An approximate, first-order expression for the 
muon intensity for > 100 GeV is [7] 



and 
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where Z^n and ^.ya' ^ire spectrum weighted mo- 
ments for production of pions and kaons, </>(£') = 
E dN/dE and (poiE^) is the intensity of primary 
cosmic-ray nucleons evaluated at the energy of the 
muon. The kinematic factors are A 



A 



0.67, 

0.25 (including the branching ratio for 

K ^iv), B^f, « 1.07 and Bk^ « 1.13. 

Zatsepin and Kuz'min pointed out long ago [8] 
the potential of a measurement of the atmospheric 
muon flux as a function of zenith angle for mea- 
suring the relative importance of kaon to pion pro- 
duction in hadronic interactions. The angular de- 
pendence arises from the denominators of the two 
terms in Eq. 1 and the numerical values of the 
critical energy parameters, e-jr ~ 115 GeV and 
CK ~ 850 GeV. In the TeV range, the kaon con- 
tribution is relatively more important near the ver- 
tical than at large angle. Thus the angular de- 
pendence of the muon flux is sensitive to the ra- 
tio Znk/Zntv- The angular dependence of the 
muon charge ratio provides information on the rel- 
ative importance of kaons and pions separated by 
charge. By fitting the data of Fig. 1 to the charge 
ratio calculated from Eq. 1, the MINOS group find 
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The ratio {Z]^k++Z^x-) I {'Z' n -k^ + "Zi n ■k-)'^'^'^ 
been kept fixed at its standard value. The fit uses 
the data of Refs [4, 10], which have data binned 
both in and cos Q. The preliminary measure- 
ment of the charge ratio in the shallow MINOS 
near detector [9] (not shown) is consistent with the 
L3+C data. The large value of the K'^ /K~ ratio 
reflects the importance of forward associated pro- 
duction (p ^ A if +) which is amplified in the 
spectrum-weighted moment by the steep primary 
cosmic-ray spectrum because the carries on 
average a significant fraction of the beam energy. 
Two other papers from MINOS are noteworthy 
as examples of the use of cosmic-ray muons for 
calibrating deep detectors. More than 20 mil- 
lion muons were measured (after cuts) in the MI- 
NOS far detector over a three-year period from 
August 2003 to August 2006. One analysis uses 
the shadow of the moon to determine the angu- 
lar resolution and absolute pointing of the far de- 
tector [11]. Both the resolution and the absolute 
pointing are 0.3° ± 0.05°. The moon's shadow is 
seen at the level of 4 ct. The shadow of the sun is 
also seen, but at a somewhat lower significance, in 
part because of bending of the parent cosmic rays 
in the solar magnetic field. 

Another paper [12] presents the analysis of sea- 
sonal variations of the underground muon rate ob- 
served in the MINOS far detector. The observed 
rate of muons is correlated with temperature by 
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where Tes is an average of the temperature 
weighted by the probability of meson production, 
which peaks at altitudes around 15 km for trajecto- 
ries near the vertical. This correlation is explained 
in the classic paper of Barrett et al. [13] as a conse- 
quence expansion of the atmosphere when temper- 
ature increases. The correlation is large for muons 
with E^^ > > where there is a competition be- 
tween decay and re-interaction of the parent pion. 
At lower energies, most pions decay before inter- 
acting for any temperature. The observed rate, R^, 
shows a seasonal variation of « ±2% and corre- 
sponds to a value of ar = 0.87 ± 0.03, consis- 
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tent with earlier measurements by MACRO [14] 
and AMANDA [15J. 

Atmospheric neutrinos 

Similar equations to Eq. 1 describe the flux of at- 
mospheric neutrinos at high energy. However, the 
kinematic factors differ in an important way. Be- 
cause its mass is close to that of the pion, the muon 
carries most of the energy in the tt — > fiv de- 
cay. The decay K ^ jiv is, more nearly symmet- 
ric. As a result, while the kinematic parameters for 
K ^ y axe nearly equal to those for — > /it, 
they are quite different for pions. In particular, 
^TTi/ ~ 0.088 as compared to w 0.67. As a 
consequence, the dominant contribution to neutri- 
nos with Ey > 100 GeV is from kaons, and the ef- 
fect of the large charge ratio for kaons has a strong 
influence on the atmospheric neutrino spectrum at 
high energy. The tendency is to harden the TeV 
neutrino spectrum and to increase the ratio v,x/i^ij.- 
It is now widely accepted that the deficit of at- 
mospheric muon-neutrinos with its pathlength and 
energy dependence is the result of neutrino oscil- 
lations. At this conference, there were only two 
papers concerning calculation of the flux of atmo- 
spheric neutrinos. Reference [16] estimates and 
tracks the various sources of uncertainty through 
the calculation in order to evaluate the systematic 
uncertainty in the flux of atmospheric neutrinos as 
a function of neutrino energy. The contribution of 
Honda et al. [17], focuses on evaluation of rela- 
tively small effects such as variation with solar cy- 
cle and the effect mountainous overburden above 
the detector Ref. [17] is based on a revised calcu- 
lation [18] that uses a new model of meson produc- 
tion in hadronic interactions [19] based on compar- 
ison to measurements of atmospheric muons. The 
new Honda et al. neutrino flux is now closer to the 
Bartol neutrino flux [16, 20] at high energy than 
the earlier calculation [21]. In both models now, 
the ratio Zj,x+ /ZpK- is large, consistent with the 
interpretation of the increase in the ratio 
in the TeV region discussed above. The corre- 
sponding effect here is that the ratio u/i) is large 
(~ 1.7) in the TeV range). Moreover, the TeV 
neutrino flux is relatively high because of the im- 
portance of forward associated production on the 



steep cosmic-ray spectrum. (An independent con- 
firmation of the relatively high atmospheric neu- 
trino flux in the TeV range comes from the work 
of Ref. [22], which assumes the best fit osciUa- 
tion parameters and unfolds the atmospheric neu- 
trino spectrum from the Super-K measurements.) 
The imphcations for atmospheric neutrinos of the 
MINOS measurement of the muon charge ratio 
have yet to be investigated in detail, however. It 
should be possible to use the MINOS data to re- 
duce the uncertainty in existing calculations of the 
flux of atmospheric neutrinos in the TeV region 
and above. 



Neutrino oscillations 

Super-Kamiokande has been restored to its full 
complement of over 11,000 50 cm photomultipli- 
ers plus an outer veto detector and has been oper- 
ating since July 12, 2006 as Super- Kamiokande- 
III. Reference [26] reviews the history of Super- 
K, which began operation in April 1996 and an- 
nounced the discovery of oscillations of atmo- 
spheric neutrinos in 1998 [27]. The phase of oper- 
ation up to the accident in November, 2001 is SK-I. 
The detector was repaired and operated with some 
5000 PMTs redistributed to provide uniform but 
sparser coverage. SK-II ran for three years, start- 
ing October, 2002. Preliminary results of SK-III 
are in agreement with SK-I and SK-II. 

A series of Super-K papers at this conference pre- 
sented preliminary results of the combined analy- 
sis of SK-I (1489 days) and SK-II (804 days). At- 
mospheric neutrino results, for example, were pre- 
sented [29] in the same format as the main SK-I 
paper [28] . The plots of zenith angle show a ratio 
of {ve + ^e)/ {Vfj. + Vfj) that is significantly higher 
than expectation for sub-GeV neutrinos from all di- 
rections, and a deficit of multi-GeV neutrinos from 
below (~10,000 km) but consistent with expecta- 
tion from above (~15 km). The angular distribu- 
tion of the electron neutrinos has the expected (no 
oscillation) shape. The results are fully consistent 
Vr osciUations with transi- 



with two-flavor 
tion probability 



sin^(2^23)xsin^ 



1.27 



(5) 



MuoNS & Neutrinos 



The dip at the first oscillation minimum in L/E is 
seen [30J for atmospheric neutrinos [31 J. and there 
is no evidence yet for three-flavor effects such as 
non-zero ^13. 

The MINOS group also presented their results 
for neutrino oscillations using the NuMI muon- 
neutrino beam from Fermilab [32]. The results, al- 
ready published, [33] show a deficit of muon neu- 
trinos in the far detector relative to the near detec- 
tor over a distance of 735 km that is consistent with 
the results of the Super-K atmospheric neutrino re- 
sult. The IVIINOS far detector can also measure 
-induced upward muons with charge separation. 
Although statistics are limited, they see a deficit 
of lower-energy neutrino-induced muons consis- 
tent with the Super-K oscillation parameters [34]. 
An interesting feature of neutrino-induced muons 
in a magnetized detector is that the charge ratio can 
be measured. The charge ratio is opposite to that 
for atmospheric muons because positive mesons 
(7r+ and K~^) decay to /i+ and (which produce 
while negative mesons give i?^ (which pro- 
duce /Lt+). 

During the time that Super-K II operated with half 
the density of PMTs as compared to Super-K I, 
new reconstruction algorithms were developed that 
allowed sensitivity similar to that of the original 
detector Now that the detector has been restored to 
its full complement of PMTs, the better algorithms 
make it possible to lower the energy threshold [35]. 
Super-K III currently is operating with 100 per cent 
trigger efficiency down to 5 MeV, which is in the 
transition region from matter dominated to vacuum 
osciUations for solar neutrinos. 

Astrophysical neutrinos (low energy) 

Neutrinos from SN1987A in the Large Magellanic 
Cloud are so far the only neutrinos detected [36, 
37] from outside the solar system. A network 
of several deep detectors continues to monitor the 
sky for bursts of neutrinos from nearby stellar col- 
lapses. New upper limits on the rate of stellar 
collapse in the Milky Way Galaxy based on non- 
observation of neutrino bursts are summarized in 
Table 5. 

It is also possible to search for a diffuse flux of relic 
neutrinos from past supernova explosions. The 



Table 2: Limits on supernova rates in the Milky 
Way Galaxy (events per year at 90% c.l.). The 
Super-K limit includes SMC and LMC. 







Experiment 


Exposure 


Limit 


S-K [44] 
LVD [45] 
Baksan [46] 


2589 d 
4919 d 
22 yr 


0.30 
0.17 
0.10 



spectrum of these reUc neutrinos peaks at a few 
MeV and falls quickly with increasing energy [38]. 
The process 

+P — >■ n -h e+ (6) 

with its relatively large cross section is the pre- 
ferred channel for this search [39]. The convolu- 
tion of the cross section, which increases with en- 
ergy, and the spectrum of reUc neutrinos may be 
above the expected background from atmospheric 
n in a window of energy from ~ 10 to 20 MeV [40]. 
Current limits from Super-K [41] are close to the 
signals expected from various models, as shown in 
Ref. [43] at this conference. The possibility [40] 
of adding Gadolinium to tag recoil nucleons from 
the process of Eq. 6 was mentioned at this confer- 
ence in Ref. [35]. A test of this method with a 2.4 
liter container of GdCla and a radioactive source is 
described in Ref. [42]. 

Astrophysical neutrinos (liigli energy) 

The most promising channel to use in the search 
for astrophysical neutrinos of high energy (> TeV) 
is neutrino-induced muons because the effective 
volume of the detector is ampUfied by the muon 
range. The average energy loss rate of a muon per 
X(g/cm^) of material traversed is 

dE 

^ = -a-bE, (7) 

where e = a/b 0.5 TeV is the char- 

acteristic energy above which stochastic losses 
(bremsstrahlung and nuclear interactions) begin to 
dominate the energy loss. The corresponding aver- 
age muon range is 

V-^M,min + e7 
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where i?^ o is the muon energy at production and 
Ef_i. min is the threshold energy of the detector for 
muons. The muon range can be several kilometers 
or more in water or ice. 

When only the muon is detected, however, there is 
only an average relation between the visible energy 
of the muon and that of the neutrino that produced 
it. The fraction of the neutrino energy carried by 
the muon varies from event to event, and the track 
is only partially contained. Moreover, for >> 
e there are large fluctuations in the amount of vis- 
ible energy deposited as the high-energy muon 
passes through the detector [23]. Nevertheless, be- 
cause the relation between energy deposition of the 
muon as it passes through the detector and its to- 
tal energy is well understood (as well as the re- 
lation between the energy of the muon and that 
of the neutrino that produced it), it is straightfor- 
ward to derive the parent neutrino spectrum from 
a measurement of neutrino-induced muons given 
sufficient statistics. From a Monte Carlo simula- 
tion of the detector response one has to identify 
a set of measurable quantities that depend on en- 
ergy deposition in the detector An unfolding pro- 
cedure can then be used to reconstruct the parent 
spectrum. The error analysis must account for the 
large fluctuations from event to event. A prescrip- 
tion using the photon density along the muon track 
as the energy-dependent observable for reconstruc- 
tion and unfolding the atmospheric neutrino spec- 
trum in IceCube is given in Ref. [24]. This is a 
natural choice given the physics of muon energy 
loss described by Eq. 7. 

The high-energy tail of the spectrum of atmo- 
spheric neutrinos constitutes the background for 
searches for neutrinos from astrophysical sources. 
Atmospheric neutrinos also serve as the calibration 
beam. Atmospheric neutrinos are sufficiently well- 
understood in the multi-TeV range so that success- 
ful reconstruction of their spectrum can be consid- 
ered as a prerequisite to any search for astrophys- 
ical neutrinos. As an example. Fig. 2 shows the 
atmospheric neutrino spectrum derived by an un- 
folding procedure from AMANDA data taken from 
2000-2003 [25]. 

An important remaining uncertainty in the atmo- 
spheric neutrino flux at high energy is the level 
of the contribution from prompt neutrinos. These 
are neutrinos from the decay of charmed hadrons 
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Figure 2: Unfolded spectrum of atmospheric neu- 
trinos by AMANDA-II [25] compared to calcula- 
tions ofRefs. [20,21]. 

which have a harder spectrum than neutrinos from 
decay of pions and kaons, which are suppressed at 
high energy because the parent mesons tend to in- 
teract rather than decay. The prompt contribution 
to the atmospheric lepton flux can be represented 
by adding a third term to the right hand side of 
Eq. 1 of the form 

A-cZn,c 
1 + BcEcos0/ec' 

where the subscript "C represents a charmed 
hadron, E is the lepton energy (fi, e, i/fj^ or z/e) and 
ec ^ 2 — 9 X 10^ for a range of charmed hadrons 
with significant leptonic branching ratios [49]. For 
muon neutrinos with > €k = 850 GeV the 
spectrum gradually steepens from E~^ '^ io E^^ '^ 
while the spectrum of prompt neutrinos continues 
to reflect the primary cosmic -ray spectrum until 
E'jy = ec ~ 3 X 10^ GeV. The crossover energy 
depends on the amount of charm production in 
hadronic interactions at high energy (Zjv (7), which 
is highly uncertain, particularly in the fragmenta- 
tion region. For a model with a significant contri- 
bution of intrinsic charm [50] the neutrinos from 
charm decay become the dominant component of 
atmospheric for E^ > 100 TeV [51]. For the 



MuoNs & Neutrinos 




o o o o o e 



this analysis 
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AMANDA-II 2000 atms. data (prelim.) 

Barr et al. atms. v + prompt atms. v 

Honda et al. atms. v + prompt atms. v 

Max uncertainty in atms. v 

Frejus 

MACRO 

AMANDA B-10 1997 diffuse 
AMANDA-II 2000 Cascades (all-flavor / 3) 
AMANDA B-1 1 997 UHE (all-flavor / 3)* 
Baikal 1 998 - 2002 ( all-flavor / 3 )* 
RICE 1999-2005 (all-flavor/ 3)* 
AMANDA-II 2000 unfolding (prelim.) 
AMANDA-II 2000-3 v„ limit 
W&B limit/2 (transparent sources) 
Full IceCube 1 yr 

' assumes a 1 :1 :1 flavor ratio at Earth 



Figure 3: Limit on the diffuse flux of astrophysical neutrinos from AMANDA-II [48]. 



much steeper spectrum of the crossover of the 
charm component is around 3 TeV. 

Diffuse search. High-energy astrophysical neutri- 
nos are expected to be produced by interaction of 
high-energy accelerated particles with gas or elec- 
tromagnetic radiation in or near the sources. Gen- 
eraUy the particle beams, and hence the produced 
neutrinos, are expected to have a harder energy 
spectrum than the background atmospheric neutri- 
nos. A standard benchmark for high-energy neutri- 
nos of extra-galactic origin is the Waxman-Bahcall 
limit [47], which assumes an E^'^ differential 
spectrum. The normalization of the Waxman- 
Bahcall limit for i/^ + Df^ at Earth after account- 
ing for oscillations is E'^dN^/dE < 2.2 x 10"® 
GeV cm'^s^^sr"^. To exceed this level would re- 
quire the existence of cosmic accelerators opaque 
to the particles they accelerate. The limit might 
also be relaxed to some extent at lower energy in 
the case of steeper source spectra. 

Figure 3 shows current limits on a diffuse flux 
of neutrinos with an differential energy spec- 
trum from AMANDA. The figure illustrates sev- 
eral points. The search labeled "this analysis" 
looks for neutrino-induced muons generated in the 



ice and rock below the detector. This is the energy 
region and the signature for which existing large 
detectors in deep water or ice are optimized. Such 
a search is limited to neutrinos with E,y < 1 PeV 
because the Earth absorbs neutrinos with higher 
energy. 

In the PeV energy region and above a diffuse sig- 
nal would be dominated by events near the hori- 
zon, where the target length is maximized with- 
out absorbing the neutrinos. Here one is gener- 
ally looking for events characterized by large and 
concentrated depositions of energy, either radiat- 
ing i/^-induced muons or cascades from interac- 
tions of or i^r in or near the detector Fig. 3 
includes the limit from the Baikal experiment in 
this energy range. The "all-flavor" limit is di- 
vided by three on the plot to make it comparable 
with the limits on i/^ alone. A preliminary value 
of 2.4 X 10"^ GeVcm^^s^^sr"^ was presented 
at this conference for the upper limit of all neu- 
trino flavors from AMANDA in the "UHE" energy 
range 10^ < E^ < 10^ GeV [52]. This is ap- 
proximately a factor of three below the Baikal limit 
and only sHghtly higher than the limit from up- 
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ward muons [48] in the energy range around 1 PeV 
where both Umits apply. 

The "diffuse" limit on any particular model spec- 
trum is obtain by convolving the spectrum with 
the detector response. In the case of the diffuse 
Umits for an spectrum, the plotted limit is a 
horizontal line on a plot of E'^dN/dE extending 
over the energy region that gives 90% of the sig- 
nal. In general, a separate limit must be calculated 
for each assumed model spectrum. Ref. [52] gives 
a useful table of models and sensitivities that spec- 
ifies which models are inconsistent at 90% confi- 
dence level with the present AMANDA UHE dif- 
fuse limit. Several early models of neutrino pro- 
duction in AGN are ruled out, including, for exam- 
ple Ref. [53], while others (e.g. [54, 55]) are stiU 
viable. 

In the case of searches for UHE neutrinos with op- 
tical detectors the signal would be characterized by 
a large amount of Ught in the detector. Since the 
events will be from above and from the sides, an 
important background is from bundles of muons 
generated by high-energy cosmic rays cascades in 
the atmosphere. Showing that this physical back- 
ground is well-understood (for example by com- 
paring simulations with data at various cut levels) 
is needed to demonstrate understanding of the de- 
tector response. A feature used to discriminate be- 
tween signal and background in Ref. [52] is the 
number of optical modules with multiple hits ver- 
sus single hits. Muon bundles from cosmic -ray 
cascades tend to produce less Ught per particle, and 
the source of the light is somewhat diffuse as com- 
pared to the intense and concentrated burst of light 
from a single particle with energy in the PeV range 
or higher. The signal would produce more multiple 
hits. 

Point sources. Particularly luminous and/or 
nearby sources of neutrinos should eventually 
emerge above the diffuse atmospheric background. 
Likely candidates are the subset of gamma-ray 
sources in which the gamma-rays are hadronic in 
origin, from decay of neutral pions produced in 
interactions of accelerated protons and nuclei in 
or near the sources. The kinematic relation be- 
tween 7r° 77 and 7r+ fJ-'^^ii provides a 
close connection between neutrinos and gamma- 
rays if the photons are not significantly absorbed 
in the sources. Examples of potential sources 



are Active Galactic Nuclei (AGNs), Gamma-ray 
Bursts (GRBs) and supernova remnants and ac- 
tive compact objects in our galaxy. Since it 
is still not known which gamma-ray sources are 
hadronic, identification neutrinos from sources of 
gamma rays (and/or electromagnetic radiation in 
other wavelengths) is the central goal of high en- 
ergy neutrino astronomy. 

As an example, it is interesting to consider Ukely 
sources of high-energy neutrinos in our local 
galaxy. In the Northern hemisphere, visible from 
Antarctica, the Cygnus region is of particular in- 
terest [56]. A systematic survey of the sensitiv- 
ity of a future kilometer-cube neutrino detector in 
the Mediterranean to potential galactic sources vis- 
ible from the North was given in Ref. [57]. (A 
more detailed discussion of the analysis is given 
in Ref. [59].) Table 3 summarizes results for those 
H.E.S.S. sources with spectra for which a break en- 
ergy (where the spectrum steepens) has been de- 
termined. The corresponding neutrino spectrum 
would steepen at a somewhat lower energy than 
the observed steepening of the gamma-ray spec- 
trum and about a factor of 40 lower than the energy 
at which the parent proton spectrum steepens. As- 
sumptions of the calculation are that the observed 
gamma-ray spectrum is entirely hadronic in origin, 
produced by interaction of an accelerated spectrum 
of protons with gas in or near the source and that 
there is no absorption of gamma-rays in the source. 
A ratio at production of ;y<. : v^, : Vr = 1:2:0 
is assumed with a flavor ratio at Earth of 1 : 1 : 1. 
The neutrino effective area of a km^ detector in the 
Mediterranean for the + 9^ channel is calcu- 
lated in some detail to obtain the expected number 
of events for source ( A^src) and background (-/Vatm) 
by convolution with the spectrum of the source and 
with the atmospheric neutrino spectrum. However, 
efficiencies for event selection and reconstruction 
are not accounted for. 

The results summarized in Table 3 nicely illustrate 
some important features of point source searches 
with kilometer-scale neutrino telescopes. The sec- 
ond and third columns of the table give the diame- 
ter of each gamma-ray source and the fraction of 
the time it is below the horizon. Typically, the 
source sizes exceed the resolution of H.E.S.S. and 
are comparable to or larger than the resolution of 
the neutrino telescope. If the characteristic neu- 
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Source name 


Dia(°) 


Vis 


e. (TeV) 






N 


-/Vatm 


A: RX J0852.0-4622 


2.0 


0.83 


1.19 


11 


104 


4.2 


21 


A:RXJ1713.7-3946 


1.3 


0.74 


1.25 


11 


41 


4.6 


8.2 


B: LS 5039 (INFC) 


0.1 


0.57 


1.01 


0.5 


2.5 


0.2 


0.5 


C: HESS J 1303-631 


0.3 


1.0 


0.21 


1.6 


11 


0.3 


2.1 


D: VelaX 


0.8 


0.81 


0.84 


16 


23 


10 


4.6 


D: HESSJ1825-137 


0.5 


0.57 


4.24 


8 


9.3 


3.7 


1.8 


D: Crab Nebula 


<0.1 


0.39 


1.72 


5.8 


5.2 


1.9 


1.1 



Table 3: TeV galactic 7-ray sources from the H.E.S.S. catalog [58] with corresponding neutrino rates cal- 
culated for 5 years operation of KM3NeT (1 km^ instrumented volume) [59]. 
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Figure 4: Preliminary sky map from Ref. [65] 
showing logio(p) for an unbinned point source 
search with AMANDA-II in 2005. 

trino break energy is > 1 TeV or higher, the signal 
to background ratio improves at higher energy, so 
the abihty to measure a signal related to energy will 
be important. For neutrino sources that steepen in 
the TeV region, the signal/background improves by 
about a factor of two if the threshold can be raised 
from 1 to 5 TeV. 

Expected rates are low, and signal/background is 
less than or comparable to unity depending on the 
size of the source. Techniques such as "source 
stacking" will therefore be important [60] to im- 
prove the significance. Similar conclusions about 
signal/background for such galactic source can be 
inferred from Ref. [61]. Hadronic models are dis- 
favored for several of the types of sources listed. 
Prime candidates are A: shell-type SNRs and C: 
TeV 7-ray sources with no counterparts at other 
wavelengths. Pulsar wind nebulae (D) and binary 
systems (B) are more often explained with electro- 
magnetic models, although hadronic models exist. 



Traditional searches use a bin size around the 
source optimized for the source size and point 
spread function of the detector. In Ref. [57, 59], 
for example, the bin size is 1.6 x \/(Tpgp -|- a%^. 
Unbinned likelihood procedures that improve sen- 
sitivity by using energy-dependence and time clus- 
tering (as well as direction) are discussed for 
ANTARES in Ref. [64], for AMANDA in Ref. [65] 
and for IceCube in Ref. [66]. Ref. [65] finds 
an improvement in sensitivity and discovery po- 
tential for AMANDA-II of 30% compared to the 
binned search. The ANTARES analysis [64] finds 
a greater improvement, up to a factor of two or 
more in some cases. Figure 4 illustrates the kind of 
confidence level map that results from an unbinned 
point source search. 

Variable sources. For sources known to be vari- 
able in electromagnetic radiation, e.g. in X- or 
7-radiation, the significance of a small number of 
neutrinos from the direction of that source could 
be greater if they occur at the same time as flares. 
Assessing the extra significance depends on the ex- 
tent to which the pattern of flaring is understood. 
Ref. [62] proposes a method in which directions 
to sources such as specific AGN known to be vari- 
able in electromagnetic radiation are searched for 
time clusters of neutrinos on various time scales. 
When a significant fluctuation above background 
is found, a check is made to see if the source is 
in a high state in electromagnetic radiation at the 
same time. If not (or if the EM data are not avail- 
able) a time-clustering algorithm is compared with 
a large set of Monte Carlo data samples from the 
selected set of sources to look for an excess over 
background. A similarly motivated approach [63] 
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looks for correlations on various time scales us- 
ing known properties of the atmospheric neutrino 
background as part of the analysis. 

The most straightforward approach to making use 
of variability in electromagnetic radiation from a 
potential neutrino source is to look for a correla- 
tion in the historical record between, for exam- 
ple, flares from blazars and times of neutrinos 
from the direction of the same sources. For co- 
incidence with a gamma-ray telescope with nearly 
continuous coverage of a large part of the sky (e.g. 
Milagro or Tibet) this is a good approach [67]. 
For telescopes with a limited field of view, how- 
ever (such as VERITAS, MAGIC, H.E.S.S.), the 
telescope will most likely be looking elsewhere 
when a neutrino signal occurs. If, as is likely, the 
neutrino events are not significantly above atmo- 
spheric background on their own, then no signal 
can be claimed. One way to address this asym- 
metry is to send an alert when a pre-specified con- 
dition is satisfied by the neutrino detector, which 
is continuously sensitive to the hemisphere below 
the detector. The gamma-ray telescope can then 
slew to the selected source and see if it is flar- 
ing. A test of such a neutrino-triggered "Target of 
Opportunity" (ToO) arrangement for a pre-selected 
set of sources was reported at this conference for 
AMANDA and MAGIC [68]. 

Another possibility is to define an alert as a group 
neutrinos from the same direction within a pre- 
selected time window for any direction in the sky. 
Such a possibility is described in Ref. [69] where it 
is proposed to send an alert to optical cameras that 
can quickly point to the direction defined by the 
group of neutrinos. In this way it might be possi- 
ble to discover the onset of an optical supernova or 
a GRB afterglow, which could elevate the signif- 
icance of the neutrino observation from a chance 
coincidence of several atmospheric neutrinos to an 
identified astrophysical neutrino event. 

Gamma-ray bursts. Neutrinos associated with 
gamma-ray bursts would have both a time tag 
and a location which would make the detection 
of even a small number of such neutrinos signifi- 
cant. Limits on neutrinos associated with gamma- 
ray bursts using AMANDA has been published re- 
cently [70, 71] and presented at this conference in 
Ref. [69]. The most sensitive search was in the 
z/^ -induced muon channel [70], which used 400 
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Figure 5: AMANDA limits on neutrinos from 
GRBs [69]. See text for discussion. 

hundred bursts reported by BATSE and IPN3 be- 
tween 1997 and 2003. No neutrinos were ob- 
served. Fig. 5 shows a comparison of limits with 
models. The three lines labeled "up-going muon" 
are limits for the models with the corresponding 
shapes. Thus the model of Ref. [72] is ruled out 
and the model of Ref. [73] is marginally incompat- 
ible with the limit. A 3cr upper limit is set at 1.3 
times the level predicted in the model of Waxman 
and Bahcall [74]. With IceCube the sensitivity for 
detection of neutrinos from GRBs will rapidly im- 
prove. IceCube is operating now with 22 strings 
and is expected to have 36 to 40 strings in oper- 
ation by the time GLAST turns on in 2008. As- 
suming GLAST will observe some 200 GRBs per 
year over the whole sky, it is estimated [69] that ob- 
servation of 70 bursts in the Northern hemisphere 
without associated neutrinos would be in conflict 
with the model of Ref. [74] at the a level. This 
level of sensitivity should be possible with an ex- 
posure equivalent to one year of full IceCube. 

Cosmogenic neutrinos 

There is now a growing consensus that the pri- 
mary cosmic-ray spectrum becomes steeper above 
5 X IQis eV [75, 76, 77]. This is generally at- 
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Neutrino Energy [eV] 

Figure 6: Figure from Ref. [80] showing upper 
limits from various experiments assuming : : 
J/,- = 1 : 1 : 1 at the detector See text for a 
discussion of this figure. 

tributed to the "GZK" effect [78, 79] of energy 
loss as particles interact with photons of the mi- 
crowave background radiation during propagation 
from sources at cosmological distances. A lack of 
high energy particles could, however, also be due 
to a lack of sources capable of accelerating par- 
ticles to energies of ^ 10^" eV. In any case, the 
number and spectrum of cosmogenic neutrinos is a 
key to the origin of the highest energy cosmic rays. 
Current measurements have not yet reached the 
sensitivity to detect cosmogenic neutrinos at the 
expected levels, as shown in Fig. 6 from Ref. [80]. 

Limits shown in Fig. 6 come from several types of 
detectors, which are sensitive to different ranges of 
energy and to different combinations of neutrino 
flavors. The Auger limit [80] is for the chan- 
nel for a period from January 2004 to December 
2006 that corresponds to one year of operation of 
the full detector. The limit is shown for the energy 
region that would generate 90% of the signal for 
an E~'^ differential spectrum, which overlaps well 
with the expected spectrum of cosmogenic neutri- 
nos. A similar limit from Hi-Res [81] is essentially 
at the same level as the Auger limit. The air-shower 



limits are based on searches for atypical horizontal 
showers as discussed in the next section. 

For uniformity, limits from other experiments are 
shown in Fig. 6 assuming an equal mixture of 
the three neutrino flavors at Earth. Limits from 
the optical detectors, AMANDA [82, 48] and 
Baikal [83], are at lower energy. Limits from ra- 
dio detectors include RICE [84] in the ice at the 
South Pole, ANITA-Lite [85], a bafloon-borne ra- 
dio detector looking for neutrinos interacting in 
the Antarctic ice sheet, and FORTE [86], search- 
ing for radio pulses from neutrino interactions in 
the Greenland ice mass. GLUE [87] looks for 
microwave signals of neutrino interactions in the 
Moon. 

There are several calculations of the spectrum of 
cosmogenic neutrinos, which vary depending on 
assumptions about the spectrum and cosmologi- 
cal evolution of the cosmic-ray sources. The band 
shown in Fig. 6 is a range based on calculations by 
two groups [88, 89]. The calculation of Ref. [89] 
was discussed at this conference in [90] for two 
different models of the primary cosmic -ray com- 
position. The result in the high-energy peak region 
{El, ^ 10^^ eV) is rather independent of the com- 
position. 

Neutrino detectors and techniques 

I conclude with a summary of the status of large 
neutrino detectors which have the primary aim of 
finding high energy (>TeV) astrophysical neutri- 
nos and identifying their sources. The discussion 
is organized by detection method. I do not in- 
clude here the densely instrumented detectors such 
as MINOS and Super-K, which are aimed primar- 
ily at study of neutrino oscillations and (in the case 
of Super-K) low energy neutrinos and proton de- 
cay. 

Because of oscillations, neutrinos from astrophys- 
ical sources are expected to consist of comparable 
numbers of all three neutrino flavors after propaga- 
tion from distant sources. At production (whether 
in the atmosphere or in an astrophysical source) 
the production of Vr is strongly suppressed rela- 
tive to and . For this reason, identification of 
T-neutrinos would be a signal of astrophysical neu- 
trinos. The signature of a tau neutrino interaction 
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Detector 


iNumoer or Uivis 


enclosed vol. (m'^) 


depth (m.w.e.) 


status 


Baikal (NT200+) 


230 


2 X 10« 


1100-1310 


operating 


AMANDA 


677 


1.5 X 10'^ 


1350-1850 


operating 


ANTARES 


900 


1 X 10"^ 


2050-2400 


2007/2008 


IceCube 


1320 


1.8 X IQS 


1350-2250 


2007 




4800 


iqs 


1350-2250 


2011 


KM3Net 


~ 10,000 




2300-3300 (NEMO) 


design study 








3000-4000 (NESTOR) 








km'^ 


1400-2400 (ANTARES site) 





Table 4: Parameters of existing and proposed neutrino neutrino telescopes in water and ice. 



ANTARES preliminary 2007 
no detector alignment 




1.1 nw 



TcTr = 49 m X £;,,(PeV). 



(9) 



Several papers at this conference focus on the phe- 
nomenology of T-neutrinos in the context of new 
experiments and the search for cosmogenic neutri- 
nos. 

Optical detectors Antares and IceCube build on 
the optical techniques of Baikal and AMANDA. 
Generally the water detectors have less scattering 
and therefore superior ability at track reconstruc- 
tion and angular resolution than ice, while being 
subject to higher background noise rates due to 
radioactivity and bioluminescence. Fig. 7 from 
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Figure 7: Figure from Ref. [91] showing angular 
distribution of reconstructed muons in 5 lines of 
ANTARES. 

is expected to become recognizable at high energy 
where the produced r-lepton has a measurable de- 
cay length. The decay length is 



Ref. [91] shows the preliminary angular distribu- 
tion of reconstructed muon tracks. The isotropic 
distribution of muons induced by atmospheric neu- 
trinos emerges from the background of downward 
muons already slightly above the horizon. This 
achievement reflects a combination of long scat- 
tering length and depth of the detector. At a depth 
of 2475 m, the intensity of penetrating atmospheric 
muons is almost an order of magnitude lower than 
at the top of IceCube (1450 m ice). At the time of 
the conference (July 2007), ANTARES had 5 lines 
operating in the ice. All 12 lines with a total 900 
of optical modules are now in place and full oper- 
ation of the detector is set to begin early in 2008. 
The status of various neutrino telescopes that use 
the optical Cherenkov technique is summarized in 
Table 4. 

IceCube is currently operating with 22 strings and 
1320 digital optical modules (DOMs) at depths of 
1450 - 2450 m in the ice at the South Pole [92]. 
Results reported at this conference were from data 
taken with the 9 string configuration that oper- 
ated during 2006. First observations of atmo- 
spheric neutrinos with IceCube-9 have been pub- 
lished [93]. The plan is to complete IceCube de- 
ployment over the next four austral summer sea- 
sons to its full size. The detector will operate and 
accumulate data while deployment is being com- 
pleted. IceCube includes a surface component, Ice- 
Top [94], which ran with 16 stations during 2006 
and is currently operating with 26 stations during 
2007. Each station consists of two ice Cherenkov 
tanks, each of which is viewed by two DOMs. Ice- 
Cube including IceTop constitutes a 3-dimensional 
air shower array that can study primary cosmic 
rays, including primary composition, from PeV to 
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Figure 8: Surface map of the coast and sea floor 
of Eastern Sicily [103]. The line marks the shore- 
line, and the highest peak is Mt. Etna. The red dot 
indicates the proposed NEMO location. 

EeV energy. AMANDA is now integrated into 
IceCube as a densely instrumented sub-array [95] 
There is a common event builder so that every 
reconstructed event in AMANDA and/or IceCube 
contains the information about hits in the detectors 
of both components. 

An important feature of the IceCube detection sys- 
tem is that waveforms of the signals are captured 
at 300 MHz. For events of high energy where 
many DOMs have multiple hits, use of full wave- 
form information can improve the reconstruction. 
For example Ref. [96] finds 5E/E - 0.34 and 
5i> ^ 0.6° for thoroughgoing tracks with 0.1 < 
< 30 PeV. 

KM3Net is a consortium of the three Mediter- 
ranean experiments (ANTARES [91], NEMO [97] 
and NESTOR [98]) to design and build a 
kilometer-scale neutrino telescope in the Mediter- 
ranean. Several papers [99, 100] present studies of 
the configuration and sensitivity of a gigaton de- 
tector in the Mediterranean Sea. Of particular in- 
terest is the comparison of three sites [102, 103] 
for their sensitivities to neutrinos with energies 
El, > 0.1 EeV, the energy range of cosmogenic 
neutrinos. In this case, neutrino-induced and 
/i^ tracks passing through the detector are of com- 
parable importance. The ly-r channel is enhanced 



by the possibility of interaction in nearby moun- 
tains from which the r-lepton can emerge and pass 
through the detector (see Fig. 8). The geometrical 
and kinematical situation are complicated by the 
competition among i^r interaction and regeneration 
in the Earth, by r energy loss and decay and by the 
configuration and response of the detector [101]. 
Accounting for these complications, the estimated 
rate is [102] 0.1 GZK neutrino interactions per km^ 
per year in a deep sea detector 

Because the Earth becomes opaque to neutrinos 
with E^, '--^PeV the angular distribution of neu- 
trinos from below can in principle be used to 
measure the neutrino cross section by varying the 
path length through the Earth, corresponding to 
the nadir angle of the event. Given a sufficiently 
large flux of neutrinos it is conceivable to sepa- 
rate neutrino cross section from energy spectrum. 
Ref. [102] give examples of energy /angular depen- 
dence in the energy range 10 PeV to 1 EeV. 

Giant air shower arrays as neutrino detectors. 

Auger and other arrays can be used to look for hor- 
izontal air showers initiated near the detector by 
a neutrino interaction with E^ > 1 EeV. Signa- 
ture of a neutrino is a horizontal shower observed 
at the ground with a large electromagnetic compo- 
nent and time structure Uke that of a normal (nearly 
vertical) cosmic-ray shower. In contrast, the back- 
ground of cosmic -ray induced horizontal showers 
consists mostly of muons with a sharp time struc- 
ture of the shower front. The electromagnetic com- 
ponent of a horizontal air shower is absorbed far 
from the detector because of the large slant depth. 




Figure 9: Diagram showing several trajectories of 
neutrino interactions that pass through the sensitive 
volume of Auger from below the horizon [104]. 



30th International Cosmic Ray Conference 



Horizontal showers from neutrinos are considered 
in two groups [101 J, those from above the hori- 
zon, which are dominated by charged current in- 
teractions of and those from below the hori- 
zon, which are mainly from v-r that interact in the 
surface of the Earth and produce a r-lepton that 
emerges from the ground and decays in the atmo- 
sphere in the field of view of the detector. The 
present limit from early operation of Auger is cited 
above (see Fig. 6). A detailed Monte Carlo simu- 
lation of Auger South to GZK neutrinos, account- 
ing for the nearby Andes mountains, is given in 
Ref. [104] and illustrated in Fig. 9. The contri- 
bution of Earth- skimming Vt is enhanced by the 
presence of the Andes mountains. In Fig. 9 lines 
labeled A and B would be from decay of r-leptons, 
while C most often would be from the interaction 
of a i^e. Class C events can also arrive from above 
the horizon. 

Radio detection of neutrinos. Given the low 
event rates expected from cosmogenic neutrinos 

in optical and air shower detectors, and in view 
of the importance of the measurement, efforts to 
find techniques that allow a much larger effective 
volume for detecting neutrinos in the EeV energy 
range are important. Using radio Cherenkov radia- 
tion from neutrino interactions in ice is one way to 
do this [105]. The Radio Antarctic Muon and Neu- 
trino Detector (RAMAND) at the Vostok station 
was the first effort to investigate radio detection in 
ice [106]. In a paper at this conference [107] some 
of the originators of the radio technique report on 
a hybrid Monte Carlo code (SIMEX) that allows 
fast simulations of radio Cherenkov radiation from 
neutrino interactions [107]. The RICE detector at 
the South Pole [108] is still in operation using the 
same technique. Studies of radio detection in ice at 
the South Pole continued in the 2006-2007 season 
with the deployment of test receivers and transmit- 
ters by AURA [109]. 

Another approach is to use a balloon-borne detec- 
tor looking down at the Antarctic ice sheet. Results 
from a prototype flight of ANITA are included in 
Fig. 6. The full ANITA detector flew over Antarc- 
tica for 35 days after launch on December 15, 
2006 [110]. Results have not yet been reported. 
Before the flight, the detector was placed above a 
large block of ice illuminated by an intense, pulsed 
electron beam at SLAC [111]. Measurements of 



10000 
1000 





















^AS^ 








le+19 le+20 le+21 le+22 le+23 

Neutrino Energy (eV) 



Figure 10: Aperture of the Square Kilometer Ar- 
ray for detection of neutrino interactions in the 
Moon [116]. 

radio pulses confirm the theory of Askaryan [112] 
on which interpretation of measurements in the 
field will be based. Aspects of propagation of ra- 
dio signals through the ice-snow-surface interface 
were checked by observing signals broadcast from 
a bore hole in the Ross Ice Shelf while ANITA 
was still in sight [113]. The possibihty of using 
air shower cores to caUbrate a radio array in the 
field was discussed in Ref. [114]. 

ARIANNA [115] is a proposal to use an array of 
10,000 antennas spread over 1000 km^ just below 
the snow surface to detect radio signals of neutri- 
nos that interact in the ice. Signals of downward 
events would be reflected by the ice-water inter- 
face. 

The history and prospects for using the Moon as 
the target and large radio telescopes on Earth as the 
detector is reviewed in Ref. [116]. Figure 10 shows 
the calculated apertures for various frequencies for 
the Square Kilometer Array [117]. The apertures 
for the original measurement with the Parkes tele- 
scope [118] are also shown. Techniques for pulse 
detection and event reconstruction are discussed in 
Ref. [119]. Although the effective areas achievable 
can be very large for SKA, as shown in Fig. 10, the 
energy threshold is well above the range for GZK 
neutrinos. 

Acoustic detection of neutrinos. This technique 
is being explored by measurements in the ice at 
the South Pole as another approach to achieving 
the effective volume needed to measure the spec- 
trum of cosmogenic neutrinos [120]. There is also 
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an acoustic test setup in Lake Baikal [121]. Ex- 
ploration of acoustic detection in ice is motivated 
by the lower noise rate in ice as compared to wa- 
ter, which should allow a lower energy threshold 
for neutrino detection. The goal is to instrument a 
sufficiently large volume to allow the detection of 
hundreds of GZK neutrinos per year. The authors 
suggest a hybrid approach [122] using acoustic, ra- 
dio and optical detectors to optimize the sensitivity 
and acceptance of the detector. 
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